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Abstract 

We consider the phenomenology of a naturally leptophobic Z-prime boson in the 1 
to 10 GeV mass range. The Z-prime's couplings to leptons arise only via a radiatively- 
generated kinetic mixing with the Z and photon, and hence are suppressed. We map 
out the allowed regions of the mass-coupling plane and show that such a Z-prime 
that couples to quarks with electromagnetic strength is not excluded by the current 
data. We then discuss possible signatures at bottom and charm factories. 
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1 Introduction 



In the past few years, the possibility of new leptophobic gauge bosons has been explored 
as a means of explaining apparent discrepancies in electroweak observables measured with 
high precision at LEP |[[ |2[], as well as an apparent high Et excess in the inclusive dijet 
spectrum at the Tevatron ||. While for the most part these anomalies have since gone 
away, the possibility remains that a Z-prime boson (Z') coupling mostly to quarks and with 
a mass smaller than mz could exist while evading experimental detection |5], ||. Given 
the assumptions that (1) the leptons are not charged under the new U(l) gauge interaction, 
and (2) the couplings to quarks are generation independent (to avoid large flavor-changing 
neutral current effects) then the normalization of the U(l) can be chosen so that the Z' 
couples precisely to baryon number. Anomaly cancellation can be achieved at the expense of 
introducing new exotic states. Two explicit examples of viable, anomaly-free models were 
presented in Refs. [f|, U, and these models presumably don't exhaust the possible ways 
in which anomalies can be cancelled. Therefore, we will set model-building issues aside 
and focus instead on the phenomenology of the Z' . This is of timely interest given the 
recent stringy suggestion that the Planck scale and weak scale might be identified |7], |j. 
In these scenarios, the dimension-5 baryon and lepton number violating operators that 
arise generically at the string scale would only be suppressed by a TeV, and hence would 
be phenomenologically lethal. Barring a higher-dimensional solution to the proton decay 
problem ||, additional gauge symmetries could provide a more prosaic, though equally 
effective, resolution. 

In Ref. ||, a specific mechanism was proposed for maintaining leptophobia in models 
with gauged baryon number, and we will adopt this mechanism here. The reason that 
leptophobia is not automatic is that the baryon number and hypercharge gauge fields mix 
via their kinetic terms 



We assume there are no Higgs fields that carry both baryon number and electroweak quan- 
tum numbers, so that mass mixing terms are not present. Below the electroweak symmetry 
breaking scale, there are separate kinetic mixing parameters for the photon and Z, which 
we will call c 7 and c^, respectively. In order that leptophobia be preserved, c 7 and cz 
must be sufficiently small at experimental energies. This can be arranged if the parame- 
ter c is forced to zero at some high scale A, so that c 7 and cz are only generated at the 
one-loop level, via renormalization group running. The boundry condidion c(A) = can 
be achieved, for example, by embedding U(1)b into a non-Abelian group, as was shown 
explicitly in Ref. |J. Here we will be more general and not assume the specific mechanism 
for achieving this boundary condition. Thus, the boundary condition, together with as- 
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sumptions (1) and (2) given above, define a class of models that we will consider further 
in the present analysis. 

In Ref. 0, the Z' mass range m T < m B < m z was studied, primarily because the 
coupling as could be taken as large as ~ 0.1 at points within this interval, without con- 
flicting with the experimental bounds. Possible high energy collider signatures were then 
considered. Here we will focus instead on Z' masses between ~ 1 GeV and my, with the 
initial goal of determining how tightly we can bound the parameter space of the model. 
Although the coupling a B cannot be as large as 0.1 within this mass range, we will show 
that current experiment does allow it to be comparable to oiem ~ 1/137. Given this result, 
we consider the possibility of detecting the Z' at charm and bottom meson factories via 
the decays of various quarkonium states which would be plentifully produced. We will not 
consider smaller values of tub, but instead refer the interested reader to the discussion in 
Ref gj. 

This letter is organized in two parts. We will first discuss the current bounds on 
the parameter space of the model. With the boundary condition on the kinetic mixing 
terms described above, both the hadronic and leptonic signatures of the Z' are completely 
determined by its mass, m B , and gauge coupling g B = ^Aho>b- Therefore, these bounds 
can be translated into boundaries of excluded regions on a two-dimensional mass-coupling 
plane. We will then consider possible discovery signals for a Z' living within these allowed 
regions. 

2 Parameter Space 

Most of the important phenomenological bounds follow directly from the Z-prime's gauge 
coupling to quarks. In addition, we take into account the small kinetic mixing effects 
by treating the mixing term in Cki n as a perturbative interaction. The Feynman rules 
corresponding to the Z — 7 and Z — Z vertices are 

-i Cl cos9 w (p 2 g^ -p^p"), (2.2) 

and 

ic z sme w {p 2 g^ -p^p v ), (2.3) 

respectively, where c 7 = cz = c above the electroweak scale, and where c = at some 
ultraviolet cutoff A. We will initially set A = mt op ~ 180 GeV, since this is probably the 
lowest scale at which the new physics responsible for the boundary condition c(A) = might 
itself remain undetected. We will describe how our results change with different choices for 
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A as needed. Note that choosing a somewhat higher value for A, for example 500 GeV, has 
only a small effect on the mixing since the dependence on A is only logarithmic. 

At any desired renormalization scale /x, we may rewrite c 7 (/i) and c^(/i) as an explicit 
function of a B by solving the one- loop renormalization group equations. These equations 
follow from the one quark-loop diagrams that connects the Z' to the 7 and Z, respectively 
i: 



8 ( \ 



2 ^Ja B a 



9n c„ 



[2N U - N d 



and 
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1 y/a B a 
18vr s 2 w c w 



[3(JV d - N u ) + 4(2iV u - N d )s\ 



(2.4) 



(2.5) 



Here c w (s w ) represents the cosine (sine) of the weak mixing angle, a is the electromagnetic 
fine structure constant, and N u (Nj) is the numbers of charge 2/3 (—1/3) quarks that are 
lighter than the renormalization scale. It is straightforward to show, for example 



c 7 (m b ) 



c 7 (m c 



0.033v^b 



cz(m b ) 



c z [m c 



0.116^05 
0.130v^ 



(2.6) 



We will use expressions like these to translate bounds on leptonic processes to exclusion 
regions on the m B -a B plane. 

The experimental bounds on the model from hadronic decays are summarized in Fig. [l|. 
Beginning with the T(1S), the new contribution to the hadronic decay width is given by 



4 
3 



a B 



a m B — m\ 



+ 



m 



T 



a m 2 B — ?7ix / 



(2.7) 



where Rf = T(T — > hadrons)/r(T — > /i + /i _ ), and the interference with s-channel photon 
exchange is included. The most stringent bound on this quantity follows from an ARGUS 
limit on the non-electromagnetic (NE) contribution to the T(1S) — > 2 jets branching 
fraction [I(J, 

BF(T(1S) -> 2 jets, NE) < 0.053 (95% CL) , 

which we find corresponds to AR^ < 2.48. This bound is stronger than the one obtained 
from the T(IS') hadronic width, discussed in Ref. |4j]. Note that we have chosen to restrict 
Fig. [I] to values of the coupling a B > 10 -3 , where direct experimental detection of the 
Z' via rare decays might be feasible. With this choice, finite width effects omitted from 
Eq. ( p.7| ) have a negligible effect on the segments of the exclusion curves shown. 

We may place additional bounds on the parameter space of the model by considering the 
hadronic decay widths of the T(2S') and T(3S") respectively. Since no direct experimental 
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Figure 1: Bounds from hadronic decays. 

bounds exist on the non-electromagnetic, two jet branching fraction, we compare Ry(2S) 
and Rr(3S) to the perturbative QCD prediction [fLlp , 



10(tt 2 -9)a£ / a s ( 3 n r , . 2a . \ 

i? = —4 ""J 1 + — - 18 - 2 + 1.161 + In 2.S 

97r a V 7r L 2 my ' ' 



where (3q — 11 — 2n//3 = 23/3. We evaluate this expression using a s (mh) as determined 



from the world average value a s (mz) = 0.119 ± 0.002 | 12| . We extract the experimental 
values of R from branching fraction data in the 1998 Review of Particle Properties ||12|| . 
This is straightforward, except in the case of the T(35*), where the branching fraction to 
ji + jjT has not been measured. We assume in this case that T(fi + fi~) is approximately equal 
to r(e + e~), which has been measured. Taking into account experimental uncertainties, we 
findAi? < 92andA# < 33 for the T(2S) and T(35) respectively, at the 95% confidence 
level. Although these bounds are weak, they nonetheless exclude some additional region of 
the parameter space immediately around the resonance masses. 

Similar bounds may be determined from the hadronic decay widths of the J/ip and 
the ijj(2S). Here, however, it is not so straigtforward to determine the standard model 
expectation. The perturbative QCD prediction for the gluonic decay width in Eq. (|2.8|) is 
derived in a nonrelativistic bound state approximation, and is therefore subject to 0(t> 2 /c 2 ) 
corrections, which are expected to be significant. Therefore, we will use the results of a 
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recent relativistic potential model analysis |T3[ as as our standard model expectation. In 
Ref. JT3"|| , the J/ip hadronic decay width was used to extract a s (m c ), yielding 0.29 ± 0.02. 
Comparing to the world average value, we find that the difference Aa s {m c ) < 0.068 can 
be tolerated, allowing two-standard deviation uncertainties. Thus, any new contribution 
to is bounded by AR < 3(Aa s /a s )R ~ 34, yielding the contour shown in Fig. |I[ We 
detemine the gluonic contribution to Rm2S) from branching fraction data in the Review of 
Particle Physics JT2j, and obtain R g = 123.6 ± 27.3. Since this is so large and uncertain, 
the bound on the model's parameter space will clearly be weak. Thus, we simply compare 
Rip(2S) to the perturbative QCD prediction, including a theoretical uncertainty comparable 
in size to the relativistic corrections in the J/ip case; we find AR < 162, yielding the curve 
shown. 

Finally, Fig. [L] displays the bound from the hadronic decay width of the Z, labelled Rz, 
which we find provides the strongest constraint from the Z-pole observables. This result 
includes the contributions from (i) direct Z' production Z — > qqZ', (ii) the Zqq vertex 
correction, and (iii) the Z — Z' mixing. These contributions were discussed in detail in 
Refs. [||, ||, using old LEP data, and here we simply include an updated bound. We will 
say nothing further on this point, since the corresponding exclusion curve is superceded by 
the others shown in Fig. |1|. 

Other bounds on the parameter depend more crucially on the kinetic mixing. We 
consider (i) the e + e~ cross section to hadrons, (ii) deep inelastic scattering, and (iii) the 
muon anomalous magnetic moment. In each case, however, we find that the constraints on 
the model are always weaker than those presented in Fig. [I]. Let us briefly consider these 
topics in turn: 

The contribution of the Z' to R = a(e + e~ — > hadrons)/<j(e + e~ — > ^ + ^~) was considered 
in Ref. ||, and was bounded by the two-standard deviation uncertainty in the experimental 
data, using a compilation of the experimental data points. While this is a reasonable 
approximation, it does not take into account that a tighter bound on any new positive 
contribution to R from a resonance effect is obtained when the central value of a given 
data point lies below the standard model prediction. Here we will take this into account, 
using the most precise measurements of R in the 5-10 GeV range obtained by the Crystal 
Ball experiment Hl4 |. Given the standard model prediction for R in the continuum region 
between the J/xjj and the T, 

R= — (l + a,/7r) « 3.54 (2.9) 
3 

we evaluate the upper bound on the difference between theory and experiment taking into 
account two standard deviation uncertainties. The tightest bound we obtained from any 
data point was AR/R < 0.05, from the measurement R = 3.31 ± 0.10 ± 0.03 ± 0.17 at 
\/s = 6.25 GeV |TJ]]. The first two experimental errors are statistical and systematic errors 
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for the given datum, while the third is an overall systematic uncertainty of 5.2%, which 
takes into account any average offset of the data. Note that within the allowed parameter 
space of Fig. [I], is not much larger than 10 -2 , and hence the Z' width is typically of 
order 10 MeV, or smaller. On the other hand, the experimental resolution at Crystal Ball 



is <Je/E = (2.7 ± 0.2)%/ JE/GeV for electromagnetically showering particles [n|, so that 
the resolution in the Z' invariant mass is comparable or larger to the Z' width. Assuming 
that m B = 6.25 GeV and as ~ 0.01 (the largest value allowed for this mass in Fig. [I]), 
we compute the contribution to AR/R by integrating the resonant and background cross- 
sections over an energy bin equal to the detector resolution, which we set equal to the Z' 
width, T = 4a B m B /9 « 28 MeV. We find 

AR/R « 0.03 

which is below the experimental bound. Since the other experimental data points present 
weaker bounds on AR/R than the one just considered, we conclude that R does not allow 
us to exclude any additional parameter space in Fig. 1. Note that at lower values of y/s 
above the charm threshold, R is not as precisely measured, and no useful bounds on the 
model can be determined. 

Deep inelastic vN scattering, parity violating eiV scattering, and the muon g — 2 provide 
only weak bounds the Z' coupling. Using the results of Ref. |J, together with the boundary 
condition described earlier, we find the corresponding exclusion regions are given by 

a B < 0.33(1 + [m B /AA7 GeV] 2 ) v N scattering (2.10) 
a B < 0.35(1 + [m B /4.47] 2 ) parity-violating e N scattering (2-H) 
a B < 1.13(m B /l GeV) 2 muon g - 2 (2.12) 

which are not even visible in Fig. [I]. Finally, we point out that resonant Bhabha scatter- 
ing places no additional bounds on the model since the nonstandard contribution to the 
amplitude is proportional to c 2 , and hence the number of events near the resonance are 
suppressed relative to the electromagnetic background by a factor of c£ ~ 10 
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Finally, we can ask how our conclusions change if the cutoff scale A is pushed to its 
largest possible value. We may use the accurate measurement the Z-hadronic width to first 
bound the mixing parameter cz(m z ); we find for a B = 0.01 that czijn z ) < 0.02. This 
corresponds to the bound A < 68 TeV. We may obtain similar bounds from consideration of 
R; however these are strongly dependent on the value of m B , as well as on the assumptions 
made in combining uncertainties from different, and often conflicting, experiments. Setting 
A to this maximum value, we find Oyimb) ~ 0.007, a factor of 2 enhancement over the 
value obtained from Eq. Q2.6|) for the same choice of a B . Clearly, this is not significant 
enough to change our qualitative conclusion that the processes involving the kinetic mixing 
in Eqs. ( |2.10| - f2~r2"D do little to constrain the parameter space of the model. 
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3 Rare Decays 



What we gather from the preceding discussion is that Fig. |T| by itself gives a reasonable 
picture of the allowed parameter space of the model. We also learn that for < tub < 
and for tub < m^, there are regions where the Z' coupling can be comparable to a^M- 
Thus, the gauge coupling need not be so small in these models as to require a separate leap 
of faith. In this section, we will assume that ~ a, and consider whether the Z' might 
eventually be detected via rare two-body decays of charm and bottom mesons. 

Since the Z' coupling to fermions is purely vectorial, the Lagrangian is charge conjuga- 
tion invariant if the Z' is C odd. This discrete symmetry forbids the decays of either the 
J /if) or T to jZ' or Z'Z' final states. Therefore, we consider instead the possible two-body 
decays of B and D mesons, as well as the decays of the lowest-lying C even quarkonium 
states, the t) c , Xc, Vb, and Xb- 

In the first case, we know that for every B or D meson decay involving a photon in 
the final state, there is an analagous process involving the Z' . The only two-body decays 
involving a photon are the various h — > S7 exclusive modes. We estimate 



T(b -> sZ') = ob + lm| m B ^ 
T(b —>■ 57) a 2 ml m 2 



— (1 + ^)(1-^) 2 (3-13) 



where « 4.3 GeV is the bottom quark mass. While this ratio is not necessarily small, 
b — > sZ' is probably not the easiest place to look for the Z' . Unlike b — > S7 which is 
discerned experimentally by study of the photon energy spectrum, b — > sZ' yields only 
hadrons in the final states, and would be overwhelmed by the larger background from 
b — > s glue |16|]. On the other hand, the contribution to the (yet unobserved) process 



b — > se + e involves the kinetic mixing, so that for as ~ a any resonance effect in the 
e + e~ invariant mass spectrum would be suppressed relative to the QED background by 



(? ~ 10 5 . The standard model prediction for the corresponding radiative decays in the D 
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meson system yield drastically smaller branching fractions, and thus, these decays are not 
likely to aid in the Z' search. 

The situation is more promising in the case of the C-even quarkonia states. For example, 
the decay r\ c — > jZ' is allowed, with 

There is an overall suppression factor of 1/4 relative to the purely electromagnetic decay 
from the squared ratio of baryon number to electric charge of the charm quark. In this case, 
one could consider rj c — > 7X, and search for a peak in the photon momentum spectrum. 
Note that the r] c branching fraction to 7 + hadrons is dominated by the decay to 7Z 7 ; 
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the decay t] c — > 7(7, where (7 is a gluon, is forbidden by color conservation, while r\ c — > 
jgg is forbidden by charge conjugation invariance. The next possibility rj c — > jggg is 
down by ~ (a^/o; s )/(27r) 4 ~ 0.001 relative to the Z' decay due mostly to phase space 
suppression, and is therefore negligible. It is simply an experimental question of whether 
single photons from other backgrounds processes can be adequately suppressed. This at 
least seems possible given that searches of exactly this type for lighter neutral gauge bosons 
have been undertaken in ir, rj and 7/ decays |17| |. A possible scenario at an e + e~ machine 



is to sit on the i[>(2S) resonance, and look for the decay chain 

ip{2S) -> 77? c -> 77X . 

One would retain events where one photon has precisely the right energy to come from 
the desired initial two body decay of the ip(2S), and then study the energy spectrum of 
the remaining photon. Exactly the same procedure could be applied to 4>{2S) — > 7Xc — > 
77X, for the various \ c states. At a charm factory with a typical beam luminosity of 
10 34 cm _2 sec _1 [18], and taking the ip(2S) production cross section to be ~ 600 nb from 
published data fL9|, we find ~ 10 4 -yZ' events per year via Xc decays, and ~ 10 3 events per 



year via r] c decays. Here we have taken the branching fraction of the \ c and i] c states to 
7Z' to be approximately 1/4 the 77 branching fractions i.e. ~ 10~ 4 . The analogous decay 
chains of the T(2S) in the 5-quark system could be studied in the same way. However, 
compared to the charmonium case, one would expect a factor of 400 reduction in the 
event rates: the production cross section for the T(2S) ||2D| is approximately two orders of 
magnitude smaller that of the ip(2S), and the jZ' branching fraction is down by a factor 
of 4 relative to the same decay in the charmonium case, due to the smaller electric charge 
of the b quark. Hence, one might still expect ~ 25 events/year from Xb decays, but the 
(yet unobserved) r] b seems less promising. 



4 Conclusions 

In this letter we have defined a generic class of naturally leptophobic Z' models, and con- 
sidered the Z' phenomenology in the 1-10 GeV mass range, a lower range than considered 
in Ref. 0. In this mass interval, decays of various quarkonia states present additional 
bounds on the Z' coupling, but new opportunities for its discovery as well. We found that 
the experimental bound on T(15') decay to two jets is primarily responsible for defining 
the allowed parameter space of the model. Bounds from the hadronic decays of the J/i/j, 
ip(2S), T(2S), and T(3S') only limit the parameter space in the immediate vicinity of the 
resonance masses; this is a consequence of larger experimental and (in the case of the char- 
monium states) theoretical uncertainties. We find that a Z' coupling as ~ «em is allowed 
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in mass intervals above and below the charmonium threshold. This opens the possibility 
of discovering the Z' in rare two-body quarkonia decays. We've suggested that perhaps 
the most interesting place to look is in the decay chain ip(2S) — > j(rj c or \c) JjZ', as 
well as in analogous decays of the T(2S). If one photon has the right energy to indicate 
an initial two-body decay to the desired quarkonium state, one could search for a peak in 
the momentum distribution of the other photon. This could provide a stunning signal of a 
light and not so weakly-coupled Z', which, given the current experimental bounds, remains 
a viable possibility. 
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